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Dynamical and statistical properties of tracer advection are studied in a family of flows produced 
by three point-vortices of different signs. Tracer dynamics is analyzed by numerical construction of 
Poincare sections, and is found to be strongly chaotic: advection pattern in the region around the 
center of vorticity is dominated by a well developed stochastic sea, which grows as the vortex system 
initial conditions are set closer to those leading to the collapse of the vortices; at the same time, the 
islands of regular motion around vortices, known as vortex cores, shrink. An estimation of the core's 
radii from the minimum distance of vortex approach to each other is obtained. Tracer transport was 
found to be anomalous: for all of the three numerically investigated cases, the variance of the tracer 
distribution grows faster than a linear function of time, corresponding to a super-diffusive regime. 
The transport exponent varies with time decades, implying the presence of multi-fractal transport 
features. Yet, its value is never too far from 3/2, indicating some kind of universality. Statistics 
of Poincare recurrences is non-Poissonian: distributions have long power-law tails. The anomalous 
properties of tracer statistics are the result of the complex structure of the advection phase space, in 
particular, of strong stickiness on the boundaries between the regions of chaotic and regular motion. 
The role of the different phase space structures involved in this phenomenon is analyzed. Based on 
this analysis, a kinetic description is constructed, which takes into account different time and space 
scalings by using a fractional equation. 

I. INTRODUCTION 

The understanding of the advection of passive tracers is of fundamental interest for many different fields, ranging 
from a pure mathematical problem, to transport or mixing related ones. One area extensively studied in the last 
decade is the so called chaotic advectionQ-Q. This phenomenon, resulting from a chaotic nature of Lagrangian 
trajectories, enhances the mixing of tracers in laminar flows, while in the absence of chaotic advection the mixing 
relies on the much less efficient mechanism of molecular diffusion. 

Chaotic advection in geophysical flows is one of the important areas of application, where the advected quantities 
vary from the ozone in the stratosphere to various pollutants in the atmosphere and ocean, or such scalar quantities 
as temperature or salinity. The interest in the geophysical flows increases the practical significance of two-dimensional 
models and more specifically the advection in the system of vortices I n addition to the large scale geophysical 

flows, 2D decaying turbulence is another example, where the inverse cascade of energy generates coherent structures 
(vortices) which dominate the evolution of the flow [lj|-|2f|. This type of problems represents only one facet of 
the interest inherent to the advection in few- vortices systems. Another facet is related to a transport of advected 
particles. It is known from different observations and numerous models, that the transport of advected particles is 
anomalous and, in one or another way, can be linked to the Levy-type processes and their generalizations pi^-pul]. 
Although these results pertain to fairly simple flows and models, there are speculations relating the chaotic advection 
in low-dimensional flows to particle dispersion in turbulent flows (see for example a discussion in 32] ). 

The interest to the chaotic advection in a three-vortex system is special not only for the reasons mentioned above. 
The three- vortex system is integrable, and its dynamics can be described in an explicit analytical form. An addition of 
a tracer (that can be regarded as another vortex of vanishing circulation) brings the number of particles in the system 
to four, which is a minimum number, from which the point vortex chaos begins the relative simplicity of the 

system permits to study anomalous transport in considerable detail. A discussion on the importance of three-vortex 
systems can be found in |^|. 

In this article, we investigate dynamical and statistical properties of the advection in flows produced by three point 
vortices with different signs. For specific conditions on both the initial position of the vortices and their strength, 
the collapse of the three vortices to a single point is then possible. We first summarize the motivations for this 
work. It is known that for systems involving a large number of vortices, the local density of vortices fluctuates, these 



fluctuations are related to situations when few vortices are close to each other and form a cluster. The number of 
vortices included in theses cluster vary, but the typical number is around 2-3 0, HE E3> H3| i while cluster involving 4 
vortices or more are much less probable. A given cluster exists only for a definite finite time At. During its existence, 
the notion of space-time structure of the cluster (or its configuration) can be introduced, and hence the influence of 
this structure on typical tracers motion and transport can be studied. Among different clusters, long time transport 
properties will be most influenced by clusters of vortices with the largest lifetime. The typical size of a cluster can be 
approximately defined by the minimum distance obtained between two vortices. The closer the vortices are to each 
other, the stronger is their mutual interaction and the longer the cluster lives. In order to change the inter-vortex 
distances it is necessary to consider a minimum of three- vortex interactions And to bring all vortices 

close together (strong interaction) the configuration between the three vortices has to be close to to a configuration 
leading to a collapse of the vortices j^. A typical lifetime estimation of these type of vortices can be assumed to be 
linked to the period of the motion of three vortices which are evolving on a close to collapse course. It is shown in 
this paper that the closer the vortices are to a collapse configuration the larger the period of the motion is, and that 
the growth of the period is exponential with respect to the closeness to the collapse configuration (see Appendix). 
We then expect that some clusters may have arbitrary long lifetime. Having in mind to shed some light on long time 
transport properties of systems involving many vortices, we decided to consider the simpler but nevertheless crucial 
case of transport properties in 3-vortex systems for close to collapse motion; this situations also provides an extreme 
situation best suited to test a possible universal transport property for three vortex flows. A simple way to define how 
close a motion is from the exact collapse course is through the deviations of the systems initial parameters from the 
conditions necessary for collapse. For instance one condition for collapse is K — (where K is related to the angular 
momentum of the system), a close to collapse motion then necessarily satisfies \K\ = e <C 1, and e is the "distance" 
from the collapse condition. Hence a special attention should be made to the notion of "close to collapse motion", 
as in the whole paper we are referring as close for a small "distance" of the system in parameter space from the 
collapse conditions. And we insist that for "close to collapse" motions in the previous sense, no collapse of the vortices 
occurs and the minimum distance of approach between the vortices is finite. In this paper the motion of the vortices is 
chosen periodic, which retains the possibility of investigating long time transport properties of these types of flows. To 
conduct our study, we use the methodology and the results of our previous works. The dynamics of tracers is analyzed 
in a spirit similar to [^3, where the structure of the advection pattern (chaotic sea, resonant islands, stochastic layers, 
coherent cores, etc) was investigated numerically and analytically for the case of three identical vortices. Transport 
properties of the advection in that particular case (where no collapse or near-collapse vortex motion is possible) were 
found to be anomalous in When vortex circulations have different signs, their dynamics may change considerably; 
the collapse phenomenon being one of the most striking examples. The motion of vortices in the vicinity of the collapse 
(when the collapse conditions are just slightly violated) was studied in j^jj, where different routes in parameter space 
all leading to the collapse were outlined. While the kinetics of advected particles in non-collapsing three-vortex flows 
was described in j^j, the situation with the advection in the near-collapse flows remained unclear. In this article we 
describe different topological structures of the advection pattern, depending on how far is the vortex system from the 
collapse conditions. We can speculate, that the characteristics of the transport of 3-vortex system can be extended 
to the case of many-vortex flows, since the advected particle finds itself quite often in the vicinity of a cluster of 3-4 
vortices, that define the advection kinetics for a fairly long time-span. For our study, we chose three different cases, 
corresponding to a specific route to collapse: a far from collapse situation, an intermediate one, and one near the 
collapse. In these three cases the collapse of the vortices is never reached and the motion of the vortices is periodic, 
which allows the study of advection for large times and a better understanding of transport in three vortex flows, as 
we consider here an extreme case with bounded motion, namely the vicinity of collapse configuration. 

In the Section 2, we present the basic equations of the point vortex dynamics, and discuss the collapse conditions 
in a three-vortex system. These results are based on the previous studies of the point vortex collapse We 
use the notations introduced in and develop some argumentation related to our choice of approaching the initial 
configuration of the vortex system corresponding to the collapse configuration. Advection equations are introduced 
in the Section 3, where we present different tools used to investigate the dynamic properties of tracers. We focus 
on Poincare sections, topology of the phase space, and trajectories, while in Section 4 we present different statistical 
results. They include velocity distributions, distributions of displacements and its moments, distribution of the 
Poincare recurrences, etc. An important point of Section 4 is to understand better how the different regions of phase 
space influence the statistical properties of trajectories of advected particles. 

On the basis of the obtained statistical information we consider kinetics of advected particles in Section 5. Fractional 
kinetics is involved in that description, and corresponding scalings and characteristic exponents are estimated on the 
basis of the results of Section 4. Motivation for the "3/2 law" of the transport is speculated, as well as the multi- 
fractal structure of kinetics. Finally, in the Conclusion we discuss different implications of the obtained results and a 
possibility to exploit them for the analysis of the advection in multi-vortex systems. 



II. NEAR-COLLAPSE VORTEX DYNAMICS 



The evolution of a system of N point vortices can be described by a Hamiltonian system of N interacting particles 
(see for instance jH3|)- The nature of the interaction depends on the geometry of the domain occupied by the fluid. 
For the case of an unbounded plane, the system's evolution writes 

dH OH 
k,zi = -2i—, 2 l =2i——,(l = l,---,N) (1) 

ozi o(kizi) 

with the Hamiltonian 

H = --!- V* kik m hx\zi - z m \ = -^-lnA , (2) 

Z7T f — ' 47T 

zi = xi + iyi is the complex coordinate of the vortex I, ki its strength and the couple (kizi,z) are the conjugate 
variables of the Hamiltonian H, to whom a new energy parameter is associated 

A = e 4 * H =Y[\ Zl -z m \ k > k ™ , (3) 

in order to simplify subsequent formulae. 

The resulting complex velocity field v is given by the sum of the individual vortex contributions: 

1 N 1 
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When zi evolves according to QJ, v provides a solution of the two-dimensional Euler equation, describing the dynamics 
of a singular distribution of vorticity 

JV 

u{z) = Y,h&(z-z l {t)). (5) 
i=i 

in an ideal incompressible two-dimensional fluid. 

The motion equations have, besides the energy, three other conserved quantities resulting from the translational 
and rotational invariance of H: 

N N 

Q + iP = Y J hz h L 2 =J2ki\zi\ 2 . (6) 

i=i i=i 

It can be easily verified, that there are three independent first integrals in involution: H, Q 2 + P 2 and L 2 , from which 
it follows, that the motion of three vortices is always integrable. An analysis of possible regimes of the motion of three 



vortices and their classification can be found in |42|, l4f 

Among the different types of motion, there is an important special case known as vortex collapse. This motion is 
available when the sum of inverse vortex circulations (harmonic mean) is zero, 

1=1 1 

and vortex positions positions are such, that the modified angular momentum K computed in the reference frame for 
which the center of vorticity is placed at the origin, vanishes 

K = ( £ ki J l2 - + p2 ) = J2 kik ™\ z * - z rf = ( 8 ) 

Vi=l / l#m 

Note, that the conditions ® and JHJ do not specify the motion uniquely, rather, they define a range of energies, for 
which the collapse is possible (once the values of ki satisfying Q are given). When both of the conditions are satisfied, 
the solutions of (QJ are singular and all three vortices collide at the center of vorticity in a finite time. Depending 



on the orientation of the vortex triangle, the collapse time t c can be either positive or negative. The first possibility, 
t c > 0, corresponds to an actual collapse, while for t c < the vortex configuration expands without bounds. These 
two cases are exact images of each other under the time-reversal symmetry, and we refer to both of them as vortex 
collapse, keeping in mind, that for t c < the collapse singularity lies backward in time. During the motion, the 
vortex configuration stays similar to the initial one, meanwhile the area of the triangle, formed by the three vortices 
grows/decreases linearly in time. We refer the reader to (4^,|42i 46] for more details. 

When the collapse conditions (|7I8|) are satisfied only approximately, the motion has a specific "near-collapse" type 
characterized by an emergence of new scales of distances and velocities, which differ significantly from the initial ones 
[4l| . A detailed analysis of the near-collapse dynamics of three vortices was performed in for the case when two 
of the vortices have the same strength. A multitude of motion regimes was found in the vicinity of collapse; different 
regimes are distinguished by the way the collapse conditions H7I8I) are violated, i.e. whether the combinations 
and K are greater, less or equal to zero. Moreover, for some of these combinations, the motion type also depends 
on the energy A of the vortex configuration; for example, when ^l/kx < and K > there exist critical energies 
A Cl and A C2 , such that in the range A Cl > A > A C2 the motion has two branches of periodic motion, for A = A C2 
the motion is aperiodic, and for A > A C2 there is only one periodic branch left. A classification of the near-collapse 
motion regimes, including the behavior of their length and time scales as the collapse is approached, can be found in 

m. 

In order to carry out a detailed study of the advection for close to collapse situations, we restrict our consideration 
to a specific one-parameter family of near-collapsing vortex systems, defined as follows: 

1. Only one of the collapse conditions, the strength condition JJJ, is violated; the second condition JHJ is satisfied, 
i.e. K = 0. 

2. The two positive vortices are identical, by an appropriate choice of time units their strength can be put to 1; 
in other words, we fix k\ = k% = 1. The circulation of the third vortex is negative, k% = —k, (k > 0). In this 
situation, the collapse happens, when the strength of the third vortex reaches a critical value k = k c = 1/2. The 
"distance" of the system from the collapse can be measured by the amount of the strength detuning S, 

6= k c - k = 1/2- k (9) 

We will be considering only the case 5 > 0. 

3. The energy of the vortex configuration A is fixed to a constant value A = 0.9. 

The first two conditions ensure the relative vortex motion to be periodic, and the choice of energy is such, that in 
the approach to collapse (6 — ► 0), the maximum inter- vortex distance grows very slowly (no noticeable changes in the 
range of 5 considered below) while the minimum inter- vortex distance rapidly approaches zero. It gives us a convenient 
"collapse in a box" setting, where vortices initially separated by a distance of order one, are brought arbitrary close 
(controlled by 5) together, which may be of interest for studies of transport in many vortex-systems. For instance, 
the ability to bring vortices closer to each other by orders of magnitude makes these 3-vortex processes a dominant 
interaction mechanism in a rare gas of vortex patches 38] . 

Before proceeding to the discussion of the advection, we will briefly summarize the results from ^lj) pertaining to 
our case. The dynamics of a three-vortex system with two identical vortices can be mapped to a one-dimensional 
Hamiltonian system describing a motion of a particle of a unit mass and zero total energy in an effective potential, that 
depends on the strength of the third vortex k, and the constants of motion A and K (see appendix). The dynamical 
variable X of this one-dimensional system is equal to the squared distance between the two positive vortices 

X = \z 2 - Zl \ 2 , (10) 

the other two distances, can be found from the expressions for A and K. And in our special case (K = 0) the motion 
is confined between two single roots of the potential, 

Xrnin ^ A <C X raax , (1 1) 

where 

k/s , 1 \ k/S 



(12) 



which implies that X{t) is a periodic function of time, and consequently the other inter-vortex distances are too, i.e. 
the relative motion of vortices in our one-parameter family is always periodic. 



As the collapse is approached (fe — ► 1/2, 8 — ► 0), the motion tends to cover all length scales: X m i n — > 0, A ma;E — > oo, 
and its period diverges as 

T~\k- 1/2S . (13) 
o 

(see the appendix for the derivation) . The influence of this behavior on the properties of advection is investigated in 
the next section. 



III. DYNAMICS OF THE ADVECTION 



A passive particle (tracer) follows the flow according to the advection equation 

z = v(z,t) (14) 

where z(t) represent the tracer trajectory, and v is the velocity field. In the case of a point vortex system, the velocity 
field is given by Eq. £Q . The incompressibility of the flow allows to write the advection equation l|T4T) in a Hamiltonian 
form: 

i = -#, *-<£ (15) 

oz oz 



where a stream function 

1 3 

*(z,z t t) = -—J2k l hi\z-z l (t)\ (16) 



1=1 

acts as a Hamiltonian. This system in non-autonomous, since the stream function depends on time through the vortex 
coordinates zj(t). The character of this dependence (periodic or not) is important for the further analysis. Below we 
will show, that although iflfijl is quasiperiodic, it can be made periodic by an appropriate coordinate transformation, 
which means, that the advection in our system has all/2 degrees of freedom Hamiltonian dynamics. 

Indeed, as was mentioned in the previous section, the relative vortex motion is periodic, i.e. the vortex triangle 
repeats its shape after a time T. This does not imply a periodicity of the absolute motion, since the triangle is rotated 
by some angle 6 during this time, see Fig. ^ In general, 9 is incommensurate with 2ir, rendering a quasiperiodic 
time dependence of Zi(t). 

Let us consider a reference frame, rotating around the center of vorticity with an angular velocity 

Q = e/T. (17) 

In this co-rotating reference frame, vortices return to their original positions in one period of relative motion T, see 
Fig. [21 their new coordinates 

I = z e~ mt , 



are periodic functions of time. 

In the co-rotating frame the advection equation retains its Hamiltonian form with a new stream function ^ which 
acquires an extra (rotational energy) term 

* = + n 2 /2\z\ 2 . (18) 
An advantage of this new frame is that VP is time-periodic: 

V{z,z,t + T) = V(z,z,t) (19) 

and well-developed techniques for periodically forced Hamiltonian systems can be used to study its solutions. 

Note, that the one-period rotation angle 8 is defined modulo 2ir, making the choice of the co-rotating frame non- 
unique. We remove this ambiguity by requiring the negative vortex to make no full revolutions around the center of 
vorticity in the co-rotating frame (as in Fig. EJ. This particular choice of O is inconsequential for the further analysis. 
The one-period rotation angles, relative motion periods and angular velocities of the co- rotating frame are presented 
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17.53 
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«„ = e/r 


-0.39 


-0.54 


-0.78 



TABLE I: Different values of O(T) and the associated quantities, the relative period T and the resulting apparent rotation 
speed. 



in the Tabled As the collapse is approached (k — > 1/2), T grows rapidly (compare to the formula j!3|0 . and the 
vortices make more and more turns per period, an acceleration of the rotation speed is also observed. 

We start our analysis of the advection by numerically constructing Poincare sections of tracer trajectories (in the 
co- rotating frame) . A Poincare section is defined as an orbit of a period-one (Poincare) map P 

Pz = z{T, zo) = e-* e z(T, z ) (20) 

where z(t,Zo) denotes a solution z(t) with an initial condition z(t = 0) = z . Plots of Poincare sections for three 
different values, k — 0.2, k — 0.3 and k — 0.41 are shown in Fig ures 13151 Vortex and tracer trajectories were 
computed using a symplectic fifth-order Gauss-Legendre scheme J^fJ. Exact conservation of Poincare invariants by 
the symplectic scheme suppresses numerical diffusion, yielding high-resolution phase space portraits. 

The Poincare sections presented in Figures 13151 show an intricate mixture of regions with chaotic and regular tracer 
dynamics, typical for periodically forced Hamiltonian systems. All three phase portraits share common features with 
the advection patterns, found in a flow due to three identical point vortices [3?il l43l | : the stochastic sea is bounded by 
a more or less circular domain, there are a number of islands inside it, where the tracer's motion is predominantly 
regular. In particular, all three vortices are surrounded by robust near-circular islands, known as vortex cores. 
Contrary to the case of three identical vortices, where the tracer dynamics is integrable for a special value of vortex 
energy (when the vortices form a steadily rotating isosceles triangle) and has a near-integrable character in the vicinity, 
the tracer motion in the near-collapse flow family considered here is always strongly chaotic, the stochastic sea remains 
a principal element of the advection pattern for any k. The phase portraits indicate, that the degree of chaotization 
increases with the approach to collapse. For instance, in the far from collapse case k = 0.2 the islands of regular 
motion inside the stochastic sea occupy a considerable area, and as k increases, their share drastically diminishes. 

A decrease in the radii of the vortex cores is of particular interest, since they are the robust structures, which also 
appear in many-vortex systems. The upper bound on the core radii can be obtained from the minimum approach 
inter-vortex distances. The minimum distance between the two positive vortices (see 111211 1 is 

/ 1 _ u\ k/26 

at the same moment the distance between the negative vortex and one of the positive two also reaches its minimum, 
which is 




R min ( 22 ) 

At this moment the vortices are collinear. Since the sum of the core radii of two vortices cannot exceed the minimum 
distance between them, we get an upper bound for the radius of the positive vortex core, Rt ore an d for that of the 
negative one, i?J ore , in terms of R" lm , R 



nun . 
2 



Rcore - mm 1 ' 1 + 2 



n— v nmin /q Q> 

-"■core — , it 2 \ Li> ) 

1 + vfc 

where we took into account, that the minimum distance between the positive and the negative vortex is shared between 
the corresponding cores depending on relative strength, and choose to define the boundary as the point between the 
two vortices with minimum speed. The core radii, measured directly from the phase portraits in Figures E] El and 
their upper bounds, obtained from \T.\l are listed in the Table ITfl 
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1.48 


1.68 


2.11 



TABLE II: Comparison of the observed size of the vortex cores r + , r with their estimations R£ re,Rcore, given byESl the 
external "radius" of the stochastic sea Rsea, and the maximum distance reached between identical vortices R™ ax . 



The presence of islands where the motion is regular, in the stochastic sea, is known to alter the transport properties 
of a physical system. This phenomenon is known as "stickiness"; when a passive particle, traveling in the stochastic 
sea, gets close to an island, it is likely to stick to this island for a while and mimic a regular trajectory of a trapped 
particle. And since with each island, a whole hierarchy of smaller islands around islands is present, the particle can 
stick for arbitrary long times, which affects on the whole the transport properties of the system. 

In stickiness has been exhibited by measuring recurrence times of particles to a given part of phase space, and 
plotting the particles position on the map with a color accordingly to their return times. Particles, that have long 
return times, all stick to a particular island, and do not jump from one to another. Taking these facts into account, we 
use another way to visualize stickiness. Indeed, sticking particles have all long coherent time behavior, which reflects 
in quantities such as their angular speed, or intrinsic speed. We decide then to compute the average intrinsic speed 
over a definite amount of time of an ensemble of particles and record it. The measured quantity is the following 

V t {m,n) = — / Vi(t)dt, (24) 

where n is the number of periods over which the speed is averaged, m keeps track of the elapsed time, and Vi(t) is the 
instantaneous speed at time t of the particle i. We then define the distribution of such averaged velocities as 

P (V; n,m) = ^J2 6 ^ V - V ^ n ^ ' ( 25 ) 

and smooth it over an interval to obtain a continuous curve. In fact, as can be observed in Fig.^J for which the speed 
of particles is plotted versus time, we can notice that after a brief period of time the distribution seems stationary, 
meaning that p(V,n, m) is independent of m and therefore we can average its value over ro, which in practice allows 
better statistics. FigureElis already informative, as we can notice some darker stripes, which means that some special 
average velocities are favored. This fits with the picture of some particles sticking to some island for a long time (at 
least > 10T here). To obtain this data we computed the trajectories over 10000 periods for a sample of 253 particles 
and recorded every n = 10 periods. We use the stationarity property and plot the distribution p versus the speed 
V for the three different cases k = 0.2, 0.3, 0.41, these are represented respectively in the figures Fig. Fig. and 
Fig. El We notice that the dark stripes observed in Fig. correspond to peaks in the density probability. In order 
to characterize the origin of these peaks, we plotted in Fig. Fig. 111! and Fig. El the position of the particles in 
the phase space contributing to the peaks in the distribution function. As anticipated, each of the observed peak 
correspond to a specific region of the phase space around some island. Concerning the influence of collapse, we notice 
that the area occupied by the contributing particles decreases as the critical condition is approached. We may then 
speculate that the transport properties of the three different system, which we discuss in the next section, may differ 
in a substantive way. 



IV. ANOMALOUS STATISTICAL PROPERTIES OF TRACERS 



Deterministic description of the motion of a passive particle in the mixing region is impossible, since a local instability 
produces exponential divergence of trajectories, and after a short time, the position of a tracer would be completely 
unpredictable. Even the outcome of an idealized numerical experiment is non-deterministic in this situation, since a 
round-off error is creeping slowly but steadily from the smallest to the observable scale. For this reason, long-time 
behavior of tracer trajectories in the mixing region is usually studied within a probabilistic approach. 
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1.589 

T > 5 • 10 5 



TABLE III: Basic values of transport properties. The average rotation speed of the tracers and the exponent related to the 
time evolution of the second moment are given. 



In the absence of long-term correlations, a kinetic description, which uses Fokker-Plank-Kolmogorov equation and 
leads to Gaussian statistics, j^lj works fairly well in many cases. Yet, in the present case, the complex topology of the 
advection pattern, illustrated by the Poincare sections in Figures EE indicates that one should anticipate anomalous 
statistical properties of the tracers in the chaotic sea. Singular zones around KAM islands usually produce long-time 
correlations, which may result in essential changes in the particle kinetics. Although in some cases these "memory 
effects" can be accounted fo r by the modification of the diffusion coefficient in the FPK equation often their 

influence is more profound |5(f[H4jp||40j. and leads to a super-diffusive behavior with faster than linear growth of 
the particle displacement variance: 

((x-(x)) 2 )~t>* (26) 

where the transport exponent fi exceeds the Gaussian value: \i > 1. 

In this section we analyze the statistical properties of tracers in the chaotic region for three vortex flow geometries, 
introduced above: far from collapse (k = 0.2, Fig. UU, intermediate (k = 0.3, Fig. QJ, and close to collapse (k = 0.41, 
Fig. EJ. A plot of a time series of the arclength versus time s,(f) = J Vi(t)dt for a set of typical tracers trajectories 
(Fig. El reveals an intermittent character of tracer motion: random pieces of trajectory are interrupted by regular 
flights, some of which are fairly long. To remain consistent with previous work, we focus our interest on the character 
of tracer rotation, and for that matter, we define its azimuthal coordinate in the center of vorticity reference frame 

0(f) = Arg z (27) 

to be a continuous function of time, i.e. 0(f) £ (— oo, oo) keeps track of the number of revolutions performed by a tracer. 

Mean advection angle (0(f)) (here ( ) denotes ensemble average) grows linearly with time: 

(0(f)) = tut, (28) 



the values of the average rotation frequency lu for the three cases are listed in Table IIIII 
The growth of the variance, 

a 2 (f) ee ((0(f) -(0(f))) 2 ) (29) 

is faster than linear for all three cases: angular tracer diffusion is anomalous. From log-log plots of er 2 (f) versus time 
in Fig. El one may see, that in order to describe the growth of the variance with a power law 

o- 2 (f)~f^ (30) 

one has to introduce different transport exponents for different time ranges. The values of these exponents, obtained 
by linear fits of corresponding parts of the graphs of Fig. El are given in Table ITTT1 Below, the first time range (with 
the exponent im) will be referred to as "short times", and the second one {(12) "long times". 

Recently [55], two types of anomalous diffusion were distinguished by the behavior of the moments, other than 
variance. The case when the evolution of all of the moments can be described by a single self-similarity exponent v 
according to 

(\x-{x)\«)~t«» (31) 
was called "weak anomalous diffusion", whereas the case when v in I'Tlll is not constant, i.e. 



(\ X -(x)\ q ) ~f^(9) 



(32) 



was named "strong anomalous diffusion". The importance of this distinction comes from the fact, that in the weak 
case the PDF must evolve in a self-similar way: 



p(x,t) = r»f(0, £ = t- v (x-{x)) (33) 

whereas non-constant v(q) in i?>'2\) precludes such self-similarity. Note, that a self-similar PDF evolution can have a 
more general form than lj33l) . with t~ v replaced by an arbitrary decaying function of time g(t): 

P(x,t)=g(t)f(g(t)(x-(x))), (34) 

which means, that if for different time decades g(t) has different asymptotics, the self-similarity exponent v will change 
from one decade to another. This variation of v with time (in particular, differences of transport exponents for short 
and long times in Table HTHl is not related to the type (strong or weak) of anomalous diffusion. 

We have performed the measurements of a set of moments of tracer angular PDF (including non-integer values of 
q) defined as: 

M 9 (t) = {\e(t)-{0{t))\ q ) (35) 
for the three vortex geometries. Time evolution of each moment was fitted by a power law: 

M q (t) ~ (36) 

separately for short, and for long times. The results are summarized in Figures ITfil ITfil fTTl (short times) and fl~8l Il9l l20l 

(long times), where the exponents n(q) are plotted versus the moment number q. In all cases, the apparent absence 
of a single linear fit indicates the presence of strong anomalous diffusion. This property was also found in ^3] (by 
comparison of the scaling properties of the central part of tracer PDF with the behavior of the variance) in a flow 
due to three vortices of equal strength. Thus, strong anomalous diffusion is a generic property of advection in three 
vortex flows. 

Our results show, that fj,(q) is well approximated by a piecewise linear function of the form: 

= ( Vq QC (37) 

yq-c if q > q c 

where c is a constant, and q c is a crossover moment number q c = c/(l — v). In where this form was introduced, 
it was found, that it fits fairly well the numerically obtained values of /i(q) in all cases of strong anomalous diffusion, 
considered there, although a theoretical example of a system with arbitrary (concave) fj,(q) was mentioned. Note, that 
deviations from l|37JI . occurring in the crossover region q sw q c , are probably a result of finite observation time, and 
the form I|37J1 might be precise in the limit t — > oo. 

As we have mentioned, the non-constant u(q) in ll32)l is incompatible with the self-similar evolution of tracer 
distribution. Let us introduce an "almost self-similar distribution" 

= / '""/(*-"*) fOT X « vt (38) 
I if x > vt 

where the exact self-similarity i'.VM) is broken only in the time-dependence of the cut-off. This modification of the 
exact self-similarity relation (j33l) takes into account the fact, that tracer speed is bounded by a maximum speed v. 
If /(£) decays sufficiently fast at infinity (e.g. like a Gaussian), the cutoff behavior is irrelevant and the moments of 
(13811 follow 13 Q , but if /(£) has a power tail, the almost self-similar distribution 138|l will give the piecewise linear 
form l|37[l for the moments. If /(£) ~ for large £, than low moments M q with q < f3 — 1 will be determined by 
the central, self-similar part of the distribution, and high moments (q > [3 — 1) by the cutoff value, 

M (t) J tm for <?</?-! m 

M ^> \ v-V-W-V if q>f}-\ m 



which is equivalent to l(37jl with 



(l-u){fi-l) (40) 



We may conclude, that the piecewise linear dependence of the exponent fi(q) on the moment number q is a signature 
of an almost self-similar evolution of tracer distribution with a long-tailed /(C)- The constant c in ll3?jl is related to 
the self-similarity exponent v and power law decay exponent (3 of /(£) by (1401 . 



Another consequence of the intermittent character of tracer motion is an anomalous distribution of recurrences of 
the Poincare map of tracer trajectories. To define recurrences, we take a region B in the chaotic sea, and register all 
returns of a Poincare map trajectory into B. The length of a recurrence is a time interval between two successive 
returns. In a system with perfect mixing, the PDF of recurrence lengths obeys a Poissonian law, provided B is small 
enough, and decay of the long-recurrence tail of the distribution is exponential for any B. Recurrence distributions 
for tracers in all three cases (k — 0.2, k — 0.3 and k = 0.41) are shown in Figures l2ll l22l l23l The plots show, that all 
distributions have long tails, indicating, that between the returns tracers are being trapped in long nights of highly 
correlated motion. The form of the graphs suggests, that long recurrences are distributed according to a power law 

P(t) ~ r~< (41) 

The values of the exponent j(S) are: 

7(0.2) = 2.2 7(0.3) = 2.4 7 (0.41) = 3.1 . (42) 

Note, that while the collapse configuration is approached, the value of the exponent increases, which may be interpreted 
as an improvement in mixing properties of the flow. This agrees with the changes in the structure of Poincare section 
(Fig. I10I12|) : the closer to collapse we get, the bigger part of the chaotic domain is occupied by a well- mixed area, 
and the smaller is the role of the singular zones around KAM islands. 

In fact one can try to find out the influence of the different islands on transport by using the distributions illustrated 
in Figures 1*7191 Indeed, each island corresponds to a specific peak. We recompute the moments of the distribution in 
the far from collapse case k = 0.2, for the modified data set, where the trajectories, corresponding to a specific peak 
are discarded. The result is presented in Fig. El We notice that the cores do affect the transport, but their influence 
is essentially visible for the high moments, while the slow particles trapped in the outer rim, are mainly responsible 
for the low moments; we also notice that we do not observe the change in slope of the strong anomalous behavior 
anymore, and conclude that the strong anomalous feature is due to the interplay of the different structures in the 
phase plane. 



V. KINETICS OF ADVECTED PARTICLES 

In some of the previous publications (see, for example, j23Is3|-jH3|) it was clearly indicated that the properties of 
anomalous transport are sensitive to phase space topology. More specifically, if we use the fractional kinetic equation 
[Hill in the form 

dPP(e,t) _ cTP{0,t) 

dtP 8\6\" [ ' 

to describe distributions P{9, t) of rotations over angle 6, then the transport coefficient V and exponents (a, (3) depend 
on the presence of different structures such as boundaries of the domain, islands, cantori, etc. The results of Section 
4 show the stickiness of trajectories of advected particles to the boundary of the domain and to boundaries of islands. 
This phenomenon is similar to what has been observed in [39] for the same-sign vortices. Our goal of this section is 
to estimate the values of the exponents a, (3. 

Figures 1101121 demonstrate stickiness of trajectories to specific structures with a filamentation of sticky domains 
along stable/unstable manifolds. In fact, different sticky domains generate different intermittent scenarios with some 
associated values of (a,/3) [EM IH^I ■ As a result, the real kinetics is multi-fractional and can be characterized by a 
set of values of (a, 0) or, more precisely, by a spectral function of (a, (3) in the same sense as the spectral function 
for multi-fractals jSfl|-[S3|. Figures 17191 show that trajectories, sticking to different structures (islands), have different 
angular velocities (compare to peaks in Figures 17191) . Due to this, different asymptotics to the distribution function 
P(9,t) and different values of (a,/3) will appear for different time intervals. In other words, for a considered time 
interval one can expect a specific "intermediate asymptotics" for P(6A) and, correspondingly, different pairs (a,/3). 
Different classes of universality for the values (a, (3) were discussed in [^j]. Below we will apply some of these results. 

Multiplying 14311 by \9\ a and integrating it over \6\ we obtain 

(IC> ~ IP (44) 

or, in the case of self-similarity the transport exponent fi from the equation 

<|0| 2 }~^ (45) 



k 


q = 2 


q > 2 


0.2 


1.6 


1.88 


0.3 


1.4 


2.0 


0.41 


1.6 


1.84 



TABLE IV: Values of the transport exponents jj, for different moments q. 



can be estimated as 

fj, = 2(3 /a (46) 

Expression 145t should be considered with some reservations since the second and higher moment may diverge. For a 
finite time t < t ma x particles reach a distance (angular rotation) 9 < 8 max , which makes all moments finite. Typically 
all 

@max — ^maxi (^^) 

and U! m ax (maximal angular velocity) can be reached only at the boundary of the domain of chaotic motion (see Figs. 
6-8). 

Using notations l.illl . i'.VHi . and 

/i(2) = fi (48) 

we can present the results for the transport exponents \i in Table IWl They are almost the same independently of how 
far is the control parameter k from its critical value k c = 1/2 (the collapse condition). For large values of q we have 
\x close to fi = 2 which corresponds to ballistic dynamics with a » j3 s» 1. 

This result can be well understood from the stickiness of trajectories to the cores (see Figs. 1101121 in the black 
color). As it follows from distributions in Figs. 17191 the particles that stick to the cores are the fastest ones, and they 
just define the large moments values. 

The value of \i for q — 2 is defined mainly by meso-structures in the middle of Figs. 1101121 (light gray). A 
typical property of these structures is existence of islands with well resolved filamentations due to the vicinity of the 
structures to a bifurcation. The latter is evident from the sharp corners of islands, which may indicate a parabolic 
typ e si ngul ar p oint [6^ |. A corresponding effective Hamiltonian, describing dynamics near a singular point, has a form 

[HllHlEllial: 

H eS = ai (AP) 2 + a 2 AQ~Q 3 (AQ) 3 (49) 

where (P, Q) are generalized momentum and coordinate and (AP, AQ) are their corresponding deviations from the 
singular point (Pq,Qo): 

AP = P - P , AQ = Q - Q . (50) 

Particularly, it may be 

Q = e, p = e. (51) 

Depending on the coefficients cij and on the meaning of variables (P, Q), which may be different from i5lTl . one can 
describe singularity due to bifurcations for different types of dynamical modes: accelerator mode f6H f64|. "blinking 
island" mode ballistic mode j^|, etc. For all these situations, universality of the Hamiltonian 149|l permits 
estimation of the exponents (a, f3) in l(4*3j) 

A trajectory that approaches the vicinity of the singular point (or, simply, a corner of the island boundary), behaves 
intermittently and escapes the near-separatrix boundary layer. The phase volume of the escaping trajectories is 

ST = 5P5Q (52) 

where SP,5Q are values AP, AQ related to the escaping particles. From l|49JI we can estimate 

^Pmax ~ SQ 3/2 (53) 

and from (1521) , (IBUll 



ST - SQ 5 / 2 



(54) 



Escaping from the boundary layer means growth of the "radial" variable 5Q with time, i.e. for an initial time interval 
SQ ~ t, and consequently, 

ST ~ i 5/2 . (55) 

From II55H we conclude for the escape probability density to leave the boundary layer at time instant t within interval 
dt: 

ip(t) oc 1/6T ~ i~ 5/2 . (56) 
It was shown in that under special conditions the exponent 7 for the trapping time asymptotic distribution 

m ~ t-T (57) 

can be linked to fractal time dimension. Moreover, 7 is related to the kinetic equation (Ii3l) as [58] 

= 7-1- (58) 

For the considered case we have (3 = 3/2. 

For the spatial distribution of particles, the simplest situation occurs when the diffusion process has Gaussian type 
and, consequently, a = 2. In the case of the presence of hierarchical set of islands, a can be defined through scaling 
properties of the island areas. In the considered situation random walk is more or less uniform but trajectories are 
entangled near stable/unstable manifolds, i.e. in the light gray areas of Figs. 1101121 That means that a ~ 2 although 
it is not exactly 2. Finally, we arrive to: 

Ai = 2/3/a~3/2 (59) 

in correspondence to observations in the Table lTvl The value ll59|l was also discussed in j^| as one of possible universal 
values for the transport exponent fi. 

We need to comment that it is not worthwhile to try to obtain \x with a higher accuracy since a specific value of /1 
has no meaning due to multi-fractal nature of transport . It is also important that we have considered such values 
of the control parameter k for which there exists a strong filamentation. That guarantees a possibility of using Eq. 
(H9l and the following analysis. 



VI. CONCLUSION 



We have considered the dynamical and statistical properties of the passive particle advection in a family of flows, 
created by three point vortices of different signs. In all three particular cases, investigated numerically, tracer advection 
was strongly chaotic: advection patterns, visualized via Poincare sections of tracer trajectories are dominated by a 
well developed stochastic sea, occupying most of the area around the center of vorticity. With the approach of the 
vortex system to the collapse configuration, the degree of tracer chaotization increases: the stochastic sea grows, 
expanding outward and consuming some of the inner resonant islands. 

The statistics of the tracers in the chaotic region is non-Gaussian. Anomalous diffusion (faster than linear growth 
of variance) with different time and space scales was found in all three cases, as well as non-Poissonian distributions 
of Poincare recurrences (with power law decay of long recurrence probability). We did not find normal transport 
regimes, if such regime exist, they are confined to narrow windows in the parameter domain. 

Transport anomalies are caused by the phenomenon of stickiness of the chaotic trajectories to the highly structured 
boundaries of the chaotic region. In the cases considered, three important types of boundaries can be distinguished: 
external border of the chaotic sea, boundaries of the resonant islands inside the chaotic sea, and boundaries of the 
vortex cores. Each of these influences various aspects of tracer statistics, analysis of their separate contributions 
shows, that the vortex cores, that rotate with the fastest rate, determine the high moments of the tracer distribution, 
while the external boundary, being the slow est, but the most sticky, dominate the low moments. 

Vortex cores appeared in simulations US Ed, E3] i their origin and sizes were derived in for a system of three 
identical vortices; particularly it was shown that the cores are the islands of stability filled by invariant curves and 
extremely thin stochastic layers. As the control parameter k approaches the collapse value k c = 1/2, the sizes of the 
vortex cores noticeably decrease. An upper bound of the core radii, obtained from the minimum distance of vortex 
approach to each other, gives a good estimation for both positive and negative vortex core size. 

Although the transport possesses multi-fractal features, it can be successively described by a fractional kinetic 
equation with characteristic exponents a ~ 2 and (3 ~ 3/2. A corresponding moments dependence is 



(60) 



The transport can be characterized by strong intermittency which manifests itself in strong deviation from llfiOII for 
higher moments, i.e. 

(|6>| 2m ) ~ i^ m > (61) 

with /i « 3/2 for m = 1 and /i w 2m for large values of m. The latter corresponds to a strong influence of ballistic 
regime of tracer dynamics. 

We note that the value /3~3/2fora~2 has also been observed for flows generated by three identical vortices ^3] 
and since this value remains for 3-vortex flows with extreme stress (vicinity of collapse) we may reasonably speculate 
that for all periodic (bounded) three vortex flows (3 <~ 3/2. We would like to point out that the present work by 
analyzing the role played in transport by the different structures involved in the flow using various techniques, and 
by confirming a typical value of the second moment exponent should be of interest for the analysis of more realistic 
and complicated systems involving many vortices and coherent structures such as geophysical fluid dynamics. 
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EXPONENTIAL PERIOD GROWTH 

In this appendix, we recall some earlier results presented in and compute an asymptotic of the period growth 
as a function of 8 = 1/2 — k. 

It has been shown for the case of three vortices with two identical ones that the relative motion of these vortices 
can been described using an one-dimensional effective Hamiltonian 

H e ff(X,X; A, K, k) = P 2 /2 + V(X) = , (62) 

with Hamiltonian equations 

X = dH ef f jdP = P , P = -8H ef f /dX , (63) 

where X — R\ is the square of the distance between the two positive vortices, and the potential V has the following 
form 

, s UK - (1 - k)X) 2 - 4k 2 Y]\(X - K) 2 - Ak 2 Y] sUk 
V(X) ee U ^ ' > ^ 2fc2 ^ '- l - , Y = {KX) l ' k . (64) 

Let us now estimate the period of the relative motion in our case. In this paper we choose a situation with K = 0, 
using then the transformation 

U = X 2 /4k 2 Y , (65) 

the potential l(64Tl becomes, 

V(U) = J^(l-fc) 2 ^-_y (17-1) , (66) 
given the fact that H e ff = 0, the motion is confined to the negative regions of the potential which leads to 



Consequently we obtain from 116 711 the boundaries for X during the motion. We then compute the period of the 
relative motion using the effective potential 



T = 2 f X2 dX = 
J Xl ^V{X) 



In the limit 8 — > (k — » 1/2), we obtain 



V(U)~-^(U-4)(U-1), (69) 



and using the inverse transformation, 

X = (tt^V^jVP-iffl „ (A 2 !/)- 1 / 45 , (70) 
we express the period in terms of U, which leads to 

T~— ) ' =dU . (71) 

As S — > (fc — » 1/2), it is the numerator in IJZlJ, which defines the asymptotic behavior, since it is decreasing function 
of U, the dominant term is from U — 1 which leads to 



T ~ -A" 1 / 2 * 5 . (72) 



Since the condition A = 0.9 < 1, is verified (and is necessary for the collapse to happen jU), we have an exponential 
growth of the period as the vortex-collapse configuration is approached. 
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PIG. 1: Vortex trajectories for k = 0.3 in absolute reference frame. Initial vortex positions A,B (two positive vortices) and C (a 
negative one) are marked with circles. Vortex triangle A'B'C' corresponds to t = to + T/2, it is congruent to the initial AABC, 
but the two positive vortices are transposed; after another half-period of relative motion the original orientation restores. Center 
of vorticity is marked by a "+". 
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FIG. 2: Same trajectories as in Fig. ^ plotted in the co-rotating frame: vortex motion is periodic. 




FIG. 3: Poincare Map of the system for the far from collapse case. The constants of motion are K = 0, A = 0.9. Vortex 
strengths are (—0.2, 1, 1). The period of the motion is T = 10.73. 




FIG. 4: Poincare Map of the system for the intermediate case. The constants of motion are K = 0, A = 0.9. Vortex strengths 
are (—0.3, 1, 1). The period of the motion is T = 17.53. 




FIG. 5: Poincare Map of the system in the close to collapse case. The constants of motion are K = 0, A = 0.9. Vortex strengths 
are (—0.41, 1, 1). The period of the motion is T = 36.86. 
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FIG. 6: Averaged speed V over 10 periods versus time. We notice that some velocities are favored. The distribution does not 
seem to be time dependent except at the very beginning. The constants of motion are K = 0, A = 0.9. Vortex strengths are 
(-0.2, 1, 1). The period of the motion is T = 10.73. 




FIG. 7: Probability distribution of V = (s(t + 10T) — s(t))/WT, where s is the arclength of a given tracer. We notice different 
peaks. The constant of motion are K = 0, A = 0.9. Vortex strengths are (—0.2, 1, 1). The period of the motion is T = 10.726. 




FIG. 8: Probability distribution of V = {s(t+ 10T) — s(t))/lOT. We notice different peaks. The constants of motion are K = 0, 
A = 0.9. Vortex strengths are (—0.3, 1, 1). The period of the motion is T = 17.53. 




FIG. 9: Probability distribution of V = {s(t+ 10T) — s(t))/WT. We notice different peaks. The constants of motion are K = 0, 
A = 0.9. Vortex strengths are (—0.41, 1, 1). The period of the motion is T = 36.86. 




FIG. 10: Location on the Poincare section of the points with definite "special" averaged velocity (compare to peaks in Fig. 0. 
In dark gray the points corresponding to 0.595 < V < 0.605. In light and lighter gray the points corresponding respectively to 
0.81 < V < 0.83 and to 1.22 < V < 1.25. In black the points with corresponding speeds V > 1.8. The constants of motion are 
K — 0, A — 0.9. Vortex strengths are (-0.2, 1, 1). The period of the motion is T = 10.73. 
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FIG. 11: Location on the Poincare section of the points with definite "special" velocity (compare to peaks in Fig. 0. In light 
gray the points corresponding to 0.545 < V < 0.565. In black the points with corresponding speeds V > 1.5. In dark gray the 
points corresponding to V < 0.52. The constants of motion are K = 0, A = 0.9. Vortex strengths are (—0.3, 1, 1). The period 
of the motion is T = 17.53. 
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FIG. 12: Location on the Poincare section of the points with definite "special" angular velocity (compare to peaks in Fig. EJ. 
In light gray the points corresponding to 0.7488 < V < 0.7515. In black the points with corresponding speeds V > 1.5. In dark 
gray the points corresponding to V < 0.489 The constants of motion are K = 0, A = 0.9. Vortex strengths are (—0.41, 1, 1). 
The period of the motion is T = 36.86. 




FIG. 13: The deviation from the mean arclength (s(t) — Vt) versus time t is plotted for 30 particles. We notice the flights 
corresponding to a particle being in a the sticky zone around an island. The constants of motion are K = 0, A = 0.9. The run 
is over 20000 periods. The average speed is V « 0.87. Vortex strengths are (—0.2, 1, 1). The period of the motion is T = 10.73. 




FIG. 14: The second order moment for the angle distribution (((#(£) — 9(t)) 2 )) versus time for the three cases k = 0.2 (solid 
line), k = 0.3 (dashed line), and k = 0.41 (dot-dashed line). We notice a change of behavior for the large times. The constants 
of motion are K = 0, A = 0.9. 



4 8 

q 

FIG. 15: The exponent qv(q) versus the moment order q for the angle distribution ((\0(t) — 0(t)\ q ) ~ t ql/( - q ') is plotted for the 
small times (t < 3 10 4 , or t < 3000T). We notice two linear behaviors: fj,(q) = 0.79q for q < 2, and fi(q) = 0.94q + Cte for 
q > 4 Vortex strengths are (—0.2, 1, 1). The period of the motion is T = 10.7. 
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FIG. 16: The exponent versus the moment order q for the angle distribution ((\0(t) — 6(t)\ q ) ~ t^ q ') is plotted for the 
small times (t < 5 10 4 , or t < 3000T). We notice two linear behaviors: n(q) = 0.62g (q < 2), n(q) = l.OOg + Cte (q > 2). The 
constants of motion are K = 0, A = 0.9. Vortex strengths are (—0.3, 1, 1). The period of the motion is T = 17.53. 
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FIG. 17: The exponent fj,(q) versus the moment order q for the angle distribution ((\0(t) — 0(t)\ q ) ~ t^ q ') is plotted for the 
short times (t < 10 5 ). We notice two linear behaviors: fj,(q) = 0.78g (q < 2), fi(q) — 0.92q — Cte (q > 2). The constants of 
motion are K = 0, A = 0.9. Vortex strengths are (—0.41, 1, 1). The period of the motion is T = 36.85. 
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FIG. 18: The exponent fx(q) versus the moment order q for the angle distribution ({\0(t) — #(t)| 9 ) ~ t^ q >) is plotted for the 
long times (t > 10 5 ). We notice two linear behaviors: fj,(q) = 0.62q (q < 2), fj,(q) — l.Olg — Cte (q > 2). The constants of 
motion are K = 0, A = 0.9. Vortex strengths are (—0.2, 1, 1). The period of the motion is T = 10.7. 
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FIG. 19: The exponent fx(q) versus the moment order q for the angle distribution ({\6(t) — #(i)| 9 ) ~ t^ q >) is plotted for the 
long times (t > 1.5 10 5 ). We notice two linear behaviors: fj,(q) = 0.75g (q < 2), fj,(q) — 1.04q — Cte (q > 2). The constants of 
motion are K = 0, A = 0.9. Vortex strengths are (—0.3, 1, 1). The period of the motion is T = 17.53. 
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FIG. 20: The exponent n(q) versus the moment order q for the angle distribution ({\9(t) — 6(t)\ q ) ~ t^ q ') is plotted for the 
long times (10 6 > t > 5 10 5 ). We notice two linear behaviors: fx(q) = 0.77? (q < 2), (i(q) = 0.82q - Cte (q > 2). The constants 
of motion are K = 0, A = 0.9. Vortex strengths are (—0.41, 1, 1). The period of the motion is T = 36.85. 




FIG. 21: Distribution of Poincare recurrences. The constants of motion are K = 0, A = 0.9. Vortex strengths are (—0.2, 1, 1). 
The period of the motion is T = 10.7. The tail presents a power law behavior whose exponent is ~ 2.2. The simulation is 
performed over 50 000 periods, statistics are made with 1137 particles. 



7 




-1 ' 1 1 1 1 — 

10000 20000 

x 



FIG. 22: Distribution of Poincare recurrences. The constants of motion are K = 0, A = 0.9. Vortex strengths are (—0.3, 1, 1). 
The period of the motion is T = 17.53. The tail presents a power law behavior whose exponent is ~ 2.4. The simulation is 
performed over 50 000 periods, statistics are made with 1012 particles. 
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FIG. 23: Distribution of Poincare recurrences. The constants of motion are K — 0, A = 0.9. Vortex strengths are (—0.41, 1, 1). 
The period of the motion is T = 36.85. The tail presents a power law behavior whose exponent is ~ 3.1. The simulation is 
performed over 50 000 periods, statistics are made with 637 particles. 




FIG. 24: Exponents versus moment order in the case k = 0.2, the difference with FigEDor FigEDis that here we did cut some 
parts of the distribution, the star sign "*" the slow motion (near zero peak corresponding to the outer rim) has been removed, 
the circle sign "o" the fast motion (fi < — 1 corresponding to the cores) has been removed. We notice two linear behavior with 
slopes 0.78 and 1, which is respectively similar to the low and high moment behavior of Fig^lor Fig^] This allows us to 
identify the individual role in transport of the different structures in the phase space. 



